The significance of the largest scale CMB fluctuations in WMAP 
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We investigate anomalies reported in the Cosmic Microwave Background maps from the Wilkin- 
son Microwave Anisotropy Probe (WMAP) satellite on very large angular scales and discuss possible 
interpretations. Three independent anomalies involve the quadrupole and octopole: 1 ) The cosmic 
quadrupole on its own is anomalous at the l-in-20 level by being low (the cut-sky quadrupole mea- 
sured by the WMAP team is more strikingly low, apparently due to a coincidence in the orientation 
of our Galaxy of no cosmological significance); 2) The cosmic octopole on its own is anomalous at 
the l-in-20 level by being very planar; 3) The alignment between the quadrupole and octopole is 
anomalous at the l-in-60 level. Although the a priori chance of all three occurring is 1 in 24000, the 
multitude of alternative anomalies one could have looked for dilutes the significance of such a pos- 
teriori statistics. The simplest small universe model where the universe has toroidal topology with 
one small dimension of order half the horizon scale, in the direction towards Virgo, could explain 
the three items above. However, we rule this model out using two topological tests: the S-statistic 
and the matched circle test. 



I. INTRODUCTION 

It is common in science that when measurements are 
improved, old questions are answered and new ones are 
raised. In this sense, history has repeated itself with the 
spectacular new Cosmic Microwave Background (CMB) 
results from the Wilkinson Microwave Anisotropy Probe 
(WMAP) [1]. The WMAP power spectrum results [2] 
are in near-perfect agreement with the cosmological con- 
cordance model in vogue (see, e.g., [3, 4]), but have sent 
cosmologists scrambling to figure out what to make of the 
detection of a high reionization optical depth and hints 
of a running spectral index [5, 6]. The "near-perfect" 
caveat above refers to the surprisingly low amplitude of 
the observed CMB quadrupole, reflecting a lack of large- 
angle correlations significant at about the 99.9% level 
[2, 5], c.f. [7, 12]. The low quadrupole amplitude seen 
by WMAP was also observed by COBE/DMR (see, e.g., 
[8], [9] and references therein), but only now, with the 
large WMAP frequency coverage and the low detector 
noise, has it become possible to show that this low CMB 
quadrupole is not significantly affected by Galactic emis- 
sion [10-12]. Compounding the puzzle, anomalies related 
to the WMAP octopole were reported in [11]. 

During years, the suppression of large-scale power 
(sometimes referred in the literature as a cutoff), have 
been used by many authors as evidence for a "small uni- 
verse" with non-standard topology (see, e.g., [13-17]). If 
so, the conventional explanation is that such suppression 
could be caused by the fact that we live in a Universe 
with compact topology, in which power on scales exceed- 
ing the fundamental cell size is suppressed - see [18] and 
references therein for a review of this subject. If the low 
multipole behavior is caused by our universe having T 1 
[19] compact topology with one dimension small relative 
to the horizon scale, then the large scale power would 



be suppressed in this particular spatial direction as illus- 
trated in Figure 1. Another possibility is hyperspherical 
topology corresponding to a closed universe [20, 21]. 

Another class of explanations that have been proposed 
involve modifying the inflationary picture to introduce 
a cutoff in the primordial power spectrum [22-26], per- 
haps linked to the spatial curvature scale [21] or string 
physics [27]. There is also a class of models that explain 
the lower multipole values by invoking the existence of 
fluctuations in a quintessence-like scalar field which can 
introduce features on scales comparable to the present 
day Hubble radius [28]. The difficulty here is to cancel 
the large integrated Sachs- Wolfe anisotropy associated 
with a high vacuum energy density. 

The goal of this paper is to go back to the WMAP data 
and investigate these large-scale anomalies in greater 
detail, to clarify what anomalies (if any) are statisti- 
cally significant and to look for other signatures of non- 
standard topologies. The rest of this paper is organized 
as follows. In Section II, we compute the significance 
of the low quadrupole, the planar octopole and their 
alignment. In Section III A, we simulate a grid of small- 
universe models to quantify whether they fit the data by 
employing a statistic that searches for topology-induced 
symmetries in the CMB sky [19]. In Section IIIB, we 
search for matched circles in the CMB sky, an idea pro- 
posed by [29], ruling out the simplest class of small- 
universe models. Finally, we present our conclusions in 
Section IV. 
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FIG. 1: The CMB maps and their S-maps (see §111 A) . Top: the all-sky cleaned CMB map from [11] is shown on the left and its S-map 
on the right. Middle: the quadrupole map (left) and its 5-map (right). Bottom: the octopole map (left) and its 5-map (right). Note that 
all three 5-maps show dark spots in the supposed direction of suppression of its original maps, around "two o'clock" . 
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TABLE I: CMB quadrupole and octopole power. 



Measurement <5T 2 2 [/1K 2 ] 


p-valuc <5T| bjK 2 ] 


Spergel et al. model 


869.7 




855.6 


Hinshaw et al. cut sky 


123.4 


0.7% 


611.8 


ILC map all sky 


195.1 


4.8% 


1053.4 


Tegmark et al. 


201.6 


5.1% 


866.1 


Cosmic quadrupole 


194.2 


4.7% 




Dynamic quadrupole 


3.6 







II. HOW SIGNIFICANT ARE THE 
"ANOMALIES" IN THE LOWER MULTIPOLES? 

A. The low quadrupole 

The surprisingly small CMB quadrupole has intrigued 
the cosmology community ever since it was first observed 
by COBE/DMR [8]. However, it was not until the pre- 
cision measurements of WMAP [10] that it became clear 
that the anomaly could not be blamed on Galactic fore- 
ground contamination. Simulations by [5, 6] have demon- 
strated that within the context of the standard inflation- 
ary ACDM concordance model, the low large-scale power 
observed is sufficiently unlikely to warrant serious con- 
cern. The WMAP team argued that the low quadrupole 
requires a fluke at the one in 143 level and the low over- 
all large-scale power (mainly quadrupole and octopole) 
require a one in 666 fluke [5] . 

In [11], it was argued that this deficit of large-scale 
power was of lower statistical significance than previously 
claimed, and this conclusion has been strengthened by 
[7, 12]. In particular, [7] gives a thorough discussion of 
this statistical issue in the context of both frequentist 
and Bayesean analysis, concluding that the WMAP re- 
sults are in fact consistent with the concordance ACDM 
model. Here, we will further elaborate the foreground- 
related point raised in [11]. The WMAP team measured 
the angular power spectrum C'e using only the part of the 
sky outside of their Galactic cut, and [11] found that this 
cut, seemingly coincidentally, eliminated the strongest 
hot/cold spots of the quadrupole and octopole. To cor- 
rectly interpret the low a priori likelihood of the WMAP 
measurements, we need to factor it as a product of two 
different probabilities, corresponding to the following two 
questions: 

1. How unlikely is the all-sky power spectrum given a 
cosmological model like ACDM? 

2. Given the all-sky CMB map, how unlikely is it that 
the Galactic cut will eliminate such a large fraction 
of the quadrupole and octopole power? 

The key point is that the location of the Galactic cut 
is determined by the orientation of the Milky Way, i.e., 
by fluctuations that have nothing to do with the CMB 



TABLE II: Observed CMB quadrupole and octopole coefficients 
in the foreground-cleaned WMAP map [11] calculated in Galactic 
coordinates in units of fiK. 



£ m Re{a ern } Imla^} 

2 10.73 0.00 

2 1 -5.87 4.26 

2 2 -13.71 -15.15 

3 -6.52 0.00 
3 1 -9.08 0.68 
3 2 21.57 1.73 
3 3 -13.64 28.79 



quadrupole and octopole. This means that although an 
unlikely coincidence associated with the second question 
may be disturbing, it is of no cosmological significance 
and should be ignored when testing cosmological models. 

Table I summarizes measurements of the quadrupole 
(column 2) and octopole (column 4) power reported in 
[2, 5, 11]. Table II gives a breakdown of the coefficients 
ag m which produce the quadrupole and octopole power 
given in row 4 of Table I, according to 

8Tf = £(£+l)C e /2ir, 

with 

m 

C l = (2£+1)- 1 M 2 - 

e=-m 

The third column of Table I shows the probability of 
the quadrupole in our Hubble volume being as low as 
observed if the best-fit WMAP team model from [5, 6] 
is correct. For the all-sky case where £Tf has a x 2 - 
distribution with 5 degrees of freedom, the probability 
tabulated is simply 

1 - 7 [5/2, (5/2)T|/855.6 M K 2 ]/r(5/2), 

where 7 and T are the incomplete and complete Gamma 
functions, respectively. According to Question 1 , the 
low cosmic quadrupole requires about a one in 20 coinci- 
dence, and the cosmic octopole is not low at all, actually 
exceeding the theoretical prediction. The one in 666 co- 
incidence from [5] thus factors roughly into a one in 20 
cosmic coincidence (Question 1) and a one in 33 Galactic 
orientation coincidence (Question 2). 

Foreground modeling remains the dominant uncer- 
tainty underlying these arguments. Although the 
detailed tests reported in [11] suggested that the 
quadrupole and octopole contributions from the inner 
Galactic plane where unimportant (after multifrequency 
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foreground subtraction) , this issue deserves further study. 
Moreover, exotic foreground processes that are not lo- 
calized to the Galactic plane, say SZ-emission from the 
Galactic halo or the local super-cluster, would not have 
been detected in these tests. If a substantial fraction of 
the all-sky power reported in Table II turns out to be 
attributable to contamination, this will of course exacer- 
bate the problem for the ACDMs cosmological model. 

B. The quadrupole-octopole alignment 

As seen in Figure 1, the quadrupole of the foreground- 
cleaned WMAP map [11] (at middle left in the figure) and 
octopole (at bottom left) both appear rather planar, with 
most of their hot and cold spots centered on a single plane 
in the sky 1 . Moreover, the two planes appear roughly 
aligned. Can we quantify the statistical significance of 
this alignment? 

A simple way to quantify a preferred axis for arbitrary 
multipoles is to think of the CMB map as a wave function 



TABLE III: Real-valued quadrupole and octopole coefficients d( m 
(defined in Appendix A) in Galactic coordinates (left) and rotated 
into their preferred frame of equation (2) (right), in units of /iK. 
The observed STi for I = 2, 3 are the rms of these numbers times 
[£(£+l)/(2u)]V2. 



£ m Galactic Rotated 



2 


_2 


-21 43 


1 °, ^9 


2 


-1 


6.03 


-0.40 


2 





10.73 


6.72 


2 


1 


-8.30 


-0.40 


2 


2 


-19.39 


28.86 


3 


-3 


40.71 


50.58 


3 


-2 


2.45 


-1.67 


3 


-1 


0.96 


-0.50 


3 





-6.52 


-13.60 


3 


1 


-12.84 


-0.27 


3 


2 


30.50 


0.71 


3 


3 


-19.29 


-20.68 



and find the axis 8 around which the angular momentum 
dispersion 

(^|(n.L)V}^m 2 K m (n)| 2 (1) 

m 

is maximized. Here ae m (n) denotes the spherical har- 
monic coefficients of the CMB map in a rotated co- 
ordinate system with its z-axis in the the 8-direction. 
In practice, we perform the maximization by evaluat- 
ing equation (1) at all the unit vectors 8 correspond- 
ing to 3,145,728 HEALPix 2 pixel centers at resolution 
nside=512, which pinpoints the maximum to within half 
the pixel spacing, about 0.06 degrees. Table II lists the 
coefficients a tm for I = 2, 3 in Galactic coordinates. We 
then use Wigner's formula (see Appendix A) to rotate 
these coefficients into other coordinate systems, which is 
much faster than repeatedly rotating the entire sky map. 
We find the preferred axes 82 and S3 for the quadrupole 
and octopole, respectively, to be 

n 2 = (-0.1145,-0.5265,0.8424), 

8 3 = (-0.2578,-0.4207,0.8698), (2) 

i.e., both roughly in the direction of (I, b) <~ (—110°, 60°) 
in Virgo. 



Under the null hypothesis that the CMB is an isotropic 
random field, the quadrupole and octopole are statisti- 
cally independent, with the unit vectors 82 and 83 being 
independently drawn from a distribution where all direc- 
tions arc equally likely. This means that the dot prod- 
uct 82 • 83 is a uniformly distributed random variable 
on the interval [— 1,1] 3 . Equation (1) does not distin- 
guish between 8 and —8 (we find a preferred axis, not 
a preferred direction), so let us quantify the alignment 
by studying the quantity |8 2 • 83 1 which has a uniform 
distribution on the unit interval [0, 1]. Equation (2) gives 
|8 2 -8 3 | w 0.9838, corresponding to a separation of 10.3°. 
An alignment this good happens by chance only once in 
1/(1 - 0.9838) w 62. In other words, the probability that 
a random axis falls inside a circle of radius 10.3° around 
the quadrupole axis is simply the area of this circle over 
the area of the half-sphere, 1/62. 

Although [11] argued that residual foreground contam- 
ination is not likely to dominate the quadrupole or oc- 
topole, it is important to keep in mind the possibility that 
some form of residual foreground contamination might 
nonetheless contribute to the alignment of the two. Note 
that the plane orthogonal to the preferred axes given by 
equation (2) is about cos -1 0.85 w 30° away from the 
Galactic plane, which is not particularly significant (re- 
quiring only a 1 in 6 coincidence). 



1 Note that this is planarity is somewhat obscured by the Aitoff 
projection. 

2 The HEALPix package is available from http://www.eso.org/ 
science/healplx/. 



This is most easily seen as follows. In a coordinate system where 
n2 = (0,0,1) and 113 = (sin cos ip, sin sin if, cos 9), the dot 
product is n2 ■ 113 = cos which has a uniform distribution, since 
the area element dil = sin 9d9dip oc dcos0 = d(n2 -113). 
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C. The planar octopole 

We just discussed the statistical significance of two 
anomalies: an intrinsic property of the quadrupole (be- 
ing low) and the alignment between the quadrupole and 
octopole. Figure 1, however, suggests the presence of 
a third anomaly intrinsic to the octopole: being unusu- 
ally planar. In contrast, the hexadecapole and higher 
multipoles seem to exhibit a more generic behavior, as 
we expect in an isotropic random field, with no obvious 
preferred axis. How unlikely is for an octopole to be so 
planar by chance? 

According to the first statistical test we performed, us- 
ing the symmetry statistic that we will describe in detail 
in the next section, this requires a one in 20 coincidence 
in the context of a Gaussian random field. Another ob- 
vious test involves the angular momentum used in equa- 
tion (1). We define a test statistic t that measures the 
maximal percentage of the octopole power that can be 
attributed to \m\ — 3, i.e., 



t = max 



|a 3 _ 3 (n)| 2 + |a 33 (n)| ; 

|a3m(n)| 2 



(3) 



and find that for the cleaned WMAP map [11] t=94%. 
Performing a large number of Monte Carlo simulations 
corresponding to an isotropic Gaussian random field 
(where the seven real- valued coefficients a 3m are simply 
independent Gaussian random variables with zero mean 
and identical variance) , we obtain values larger than this 
about 7% of the time. Redefining t as the ratio of the 
angular momentum dispersion to the total power, i.e., 
t = (V'|(n-L) 2 |i/')/('i/>|'i/')) gives a similar significance level. 
In both of these tests, the planar nature of the octopole 
is reflected by the dominant contribution from |m| = 3 to 
the octopole rms in the last column of Table III. In con- 
trast, the quadrupole is not significantly planar according 
to these same statistical tests. 



III. TOPOLOGICAL SIGNATURES 

In the previous section, we investigated the anoma- 
lies reported in the maps from the WMAP satellite on 
very large angular scales. Although the a priori chance 
of all these three anomalies occurring is 1 in 24000, the 
multitude of alternative anomalies one could have looked 
for dilutes the significance of such a posteriori statis- 
tics. It is important to point out, however, that in 
the simplest small Universe model where the Universe 
has toroidal topology with one small dimension smaller 
than the horizon scale, one would expect all three of the 
anomalies discussed above. This means that in the con- 
text of constraining cosmic topology, these three effects 
are not merely three arbitrary ones among many. 

Two completely separate observable signatures of small 
universes have been described in the literature. The first, 



known as the ^-statistic [19, 30], assumes that modifica- 
tions of the CMB power spectrum are caused by the fact 
that only certain fluctuation modes are allowed by the 
boundary conditions, just as the fundamental tone and 
its overtones give information about the geometry inside 
a flute. The second, known as circles-in-the-sky [29], is a 
purely geometric effect in which a given space-time point 
may be observable by us in more than one direction. 

These two signatures are both independent and com- 
plementary. The information about the first comes from 
the largest angular scales in the CMB, while informa- 
tion about the second comes from smaller-scale CMB 
patterns, as well as distant objects such as quasars that 
could potentially be multiply imaged (see, e.g., [31, 32]). 
Detection of a signature of the second kind could provide 
smoking gun evidence of a small universe, being indepen- 
dent of any assumptions about cosmological perturbation 
theory and potentially providing high statistical signifi- 
cance. 



A. Searching for symmetries in the sky 

In this subsection, we investigate whether the WMAP 
data exhibits the type of large-scale symmetries predicted 
by small-universe models [19, 30]. 



1. The S -Statistic 

We perform our tests by computing the function S(hi) 
defined by [19, 30]: 



N, 



-E 



ST ST . 



(4) 



where N pix is the number of pixels in the cleaned map 
and hij denotes the reflection of hj in the plane whose 
normal is n,, i.e., 



n. 



hj - 2(hi ■ hj)hi. 



(5) 



S(h) is a measure of how much reflection symmetry there 
is in the mirror plane perpendicular to n. The more 
perfect the symmetry is, the smaller 5(n) will be. 

We refer to the map of S(hi) over the whole sky as an 
S'-map. This map is a useful visualization tool and gives 
intuitive understanding of how the S-statistic works. S- 
maps are discussed in more detail in Section III A 2. 

For a given sky map, we compute a test statistic S 
defined as the minimum value of the S'-map. As ex- 
plained and illustrated in [19], small-universe models 
with toroidal topology exhibit symmetries that give them 
on average smaller values of So than infinite universes, 
which makes the S-statistic useful for constraining the 
cosmic topology. 




FIG. 2: Cumulative histograms of S . Top-left: S-test for 
the cleaned WMAP map. Top-right: S-test for octopole alone. 
Bottom-left: S-test for the sum of the quadrupole and octopole. 
Bottom-right: S-test for hexadecapole alone. Curves show the frac- 
tion of 500 simulated maps that have So below the given value. A 
small Universe should give a small So-value, but the observed value 
of SWMAP (vertical line) is seen to be significantly smaller than 
expected in an infinite universe only for the octopole case. 



2. S-maps from the WMAP data 

The S-map of our WMAP cleaned map from [11] is 
shown in Figure 1 (top right), and shows a striking pair 
of dark spots on the same preferred axis that we identified 
in Section II B, equation (2), at "two o'clock" and its an- 
tipode. Such an isolated minimum in the S-map is char- 
acteristic of a T 1 -model, where the Universe is a 3-torus 
with only one direction small relative to the cosmic hori- 
zon. Figure 1 also shows the S-maps of the quadrupole 
(middle right) and the octopole (bottom right) compo- 
nents, and both are seen to independently show dark 
spots along this special axis. In contrast, the S-maps 
of the hexadecapole and higher multipolcs do not have 
this property. 



3. What is the probability that an infinite Universe posseses 
such symmetry? 

To quantify the statistical significance of this symme- 
try axis, we produced 500 all-sky Monte Carlo simula- 
tions with HEALpix resolution nside=16, with monopole 
and dipole removed, adopting £ max —128, and using the 
WMAP ACDM fiducial power spectra of [5]. We then 
use equation (4) to compute the S-map and S for each 
of those 500 simulations. Figure 2 (top left) shows a cu- 
mulative histogram of S for these 500 simulations. For 
comparison, the vertical line corresponds to the value of 
s wmap for the real WMAP data processed in the same 
way. The fact that S™ MAP for the real data falls near 
the middle of the distribution means that although the 
dark spots in seen in Figure 1 (top right) may appear 
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FIG. 3: Expected and observed So -values for the octopole map 
alone (top) and for the sum of the quadrupole and octopole maps 
(bottom). The solid line shows the mean of the Monte Carlo sim- 
ulations and the yellow (or grey in BW) band shows the 1 — a 
spread. 120 simulations per il^-value were performed for R x < 1 
and 30 per R x -value for R x > 1. The horizontal lines represent the 
observed values S™ MAP . 



visually striking, symmetries at this level are not at all 
unlikely to happen by chance in an infinite universe. In 
other words, the WMAP data is perfectly consistent with 
a standard infinite universe as far as our S-test is con- 
cerned. We obtain similar results when analyzing the 
octopole alone (top right), the hexadecapole alone (bot- 
tom right), various higher multipoles alone, as well as the 
sum of quadrupole and octopole (bottom left). For the 
case of the octopole (top right), however, we find the high 
degree of symmetry to be marginally significant, with an 
So-value as low as observed by WMAP occurring in less 
than 5% of the simulations. If the octopole had all its 
hot and cold spots centered on a single plane, it would 
have perfect reflection symmetry about this plane, so this 
one in 20 anomaly is simply another manifestation of the 
planar nature of the octopole discussed in Section II C. 

Note that there is nothing to be learned from applying 
the 5-test to multipoles I =0, 1 or 2 on their own, since 
(apart from pixelization effects) they all give S = 0. The 
monopole exhibits perfect symmetry around any axis, so 
its S-map is identically zero. The dipole has perfect re- 
flection symmetry in all axes perpendicular to the dipole 
direction, so its S-map equals zero on a great circle. The 
quadrupole is a quadratic function given by a symmetric 
3x3 matrix, and therefore has perfect reflection symme- 
try about the three orthogonal eigenvectors of this ma- 
trix. This is illustrated in Figure 1, where the quadrupole 
S-map (middle right) exhibits six zeroes corresponding to 
the pair of hot spots, the pair of cold spots and the pair 
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of saddle points in the quadrupole map (middle left). 

4- Making Toroidal Fake Skies 

Above we found that the S-statistic was consistent 
with an infinite universe. Since a small universe could 
nonetheless explain all three of the anomalies of Sec- 
tion II in one fell swoop, let us investigate whether small 
universes provide a better fit to the data than an infinite 
universe does or, conversely, whether the 5-statistic rules 
out some interesting class of small universe models. 

In a toroidal universe, only wave numbers that are har- 
monics of the cell size are allowed. Therefore, we have a 
discrete k spectrum [33, 34] 

Rh ( -Rx ' -Rj/ ' Rz) ^ ^ 

where p x ,p y and p z are integers, and R x , R y and R z are 
the sizes of the cell. It was shown in [19] that toroidal 
universes can be simulated by 

ST 4 
— (0, 4>) oc [9i cos(2tt7) + 32 sin(27T7)] ) 

(7) 

where 




<?i and (72 are two unit Gaussians random variables, and 
n is the spectral index 4 . We generate our simulations 
with n = 1 since this provides a good fit to the observed 
CMB power spectrum on the large angular scales that 
matter here (our tests below use only the information at 
the lowest multipoles). The symmetry patterns expected 
in T 1 universes have been shown to appear for broad 
range of n- values (actually, n <3) [19]. Since we are 
focusing on T 1 universes, we wish to set R y = R z = 00. 
In practice, we set R y — R z — 3, since this is more 
numerically convenient and essentially indistinguishable 
in practice. We add no noise to the simulations, since 
WMAP detector noise is completely negligible for I < 20. 

5. If we live in a T 1 , can we constrain the cell size R x ? 

For the tests in this section, we create two band-pass 
filtered versions of the foreground-cleaned WMAP map, 



4 A more realistic Monte Carlo can be done by adding to each 
toroidal simulation a map that has only the ISW contribution. 
A power suppression and enhanced symmetry at low multipoles 
would remain. 



retaining only £ = 3 and only 2 < I < 3, respectively. 
We then degrade the resolution of these two maps to 
HEALpix resolution nside~16. After generating our sim- 
ulated small universe maps as above, we band-pass filter 
and process them in the same way as the real data. Since 
we wish to test only for symmetry properties, we nor- 
malize the two types of Monte Carlo maps to have the 
same rms fluctuations as the two corresponding filtered 
WMAP maps. 

Fixing a cell size, we constructed a simulated CMB 
map as described above and used equation (4) to obtain 
an S'-map from which we extracted its minimum value 
S - Repeating this procedure 120 times, we obtain the 
probability distribution of S for a given cell size. When 
we repeat this procedure for different cell sizes we are 
able to construct Figure 3. 

In Figure 3, we show the probability distribution of 
S obtained from Monte Carlo simulations for cell sizes 
0.1 < R x < 3. Comparing this with the horizontal 
solid line (which represents the 5* -value extracted from 
the WMAP data, S™ MAP ) shows what range of Re- 
values WMAP favors. The octopole alone (top panel) 
is seen to favor a small Universe with R x < 0.5 at the 
1 — a level, since it has near reflection symmetry about 
the preferred plane discussed in Section II C. However, 
the quadrupolc-octopole combination (bottom panel) is 
seen to disfavor small universes, preferring R x > 1. 
We repeated the same test for all individual multipoles 

1 = 4, ...,10 and for the ranges 2 - 10, 11 - 20 and 

2 — 20, finding no significant preference for small uni- 
verses. Combining the evidence of all these tests, we 
therefore conclude that the WMAP S-test does not fa- 
vor the small universe hypothesis, preferring a cell size 
at least as large as our Hubble Volume, R x > 1. 



B. Searching for circles in the sky 

The circles method [29] involves matching circles in the 
CMB sky. The CMB Surface of Last Scattering (SLS) im- 
aged by WMAP is a sphere with us at its center. The sim- 
ple case of a toroidal Universe is mathematically equiv- 
alent to a perfectly periodic Universe, with copies of us 
on an infinite rectangular lattice. Since each copy of us 
is surrounded by its own spherical SLS, and the intersec- 
tion of two spheres a circle, it follows that if the lattice 
spacing is smaller than the diameter of the sphere, we 
can observe the same circular part of the SLS in opposite 
directions in the sky. In the case of a sufficiently small 
spacing, there will be several such matched circles, all 
concentric. As shown in [29], this result is general: any 
compact topology generates matched circles. For more 
complicated cases, these circles will typically no longer 
be opposite from one another as mirror images, but can 
differ in location, size, rotation and parity. 
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1. Real-world complications 

In the original discussion of this effect, [29] considered 
the idealized case where the entire CMB signal came from 
the SLS and would look the same from different vantage 
points. In practice, however, we need to consider three 
important departures from these assumptions: 

1. The late integrated Sachs- Wolfe (LISW) effect is 
generated not on the SLS but later on, by the grav- 
itational potential that the CMB photons traverse 
en route to us, so this contribution to the WMAP 
map will not match between paired circles. 

2. Whereas the CMB fluctuations generated by SLS 
density fluctuations and gravitational potential 
fluctuations look the same from any vantage point, 
those generated by the Doppler effect from ve- 
locity perturbations do not. Rather, they have 
a dipole structure, which means that what looks 
like a hot spot from one vantage point looks like 
a cold spot when viewed from the opposite direc- 
tion. The Doppler contribution to CMB fluctua- 
tions therefore does not match between paired cir- 
cles, and will in fact be strongly anti-correlated for 
small matched circles (whose features we view from 
nearly opposite directions). 

3. The contribution to CMB maps from foreground 
contamination and detector noise will not match 
between paired circles. 

To minimize the effect described in the third item, 
we based our analysis on the foreground-cleaned WMAP 
map from [11], which minimizes the combined fluctua- 
tions of foregrounds and noise. Although residual fore- 
ground contamination may prevent us from finding a pair 
of small matched circles hiding in the Galactic plane, the 
above-mentioned axis that we are most interested in is 
fortunately nowhere near this plane. 

We tackle the effects described in items 1 and 2 by 
filtering the WMAP map before performing the circle 
search. Since the LISW effect is important only for 
the very largest multipoles, we can high-pass filter the 
WMAP map by zeroing all multipoles with I < 8, there- 
fore removing the contributions from 1. The Doppler ef- 
fect, on the other hand, is dominant only in the troughs 
between the acoustic peaks (it is mainly this effect that 
prevents the CMB power spectrum from dropping all the 
way down to zero between the peaks). Therefore, we can 
remove the contributions from 2 by zeroing all multipoles 
with I > 300 as well. If we were to detect a suspected 
matched circle pair, a powerful subsequent consistency 
test would be whether an i = 301 — 400 band-pass filtered 
map (probing roughly the first acoustic trough) displays 
anti-correlation between the two circles. 

The issue of how to best deal with the above-mentioned 
real-world complications for circle matching is an inter- 
esting and challenging data analysis problem, worthy of 
a future paper in its own right. Our analysis here should 
simply be viewed as a first attack on the problem, and 



the forthcoming analysis by [37] , in preparation, will un- 
doubtedly improve it. One challenge is that band-pass 
filtering is spatially non-local, potentially reducing de- 
tectability by smearing out matched circles and by su- 
perimposing signal from neighboring sky regions on the 
circles. A second challenge is how to optimally deal with 
the fact that the WMAP map is pixelized. This means 
that there generally will not be a pixel centered exactly at 
the desired reflected position in equation (10) below 
in this paper, we simply use the closest pixel. This 
pixelization problem is exacerbated by a desire to per- 
form the circle search at as low resolution as possible to 
avoid the computations becoming prohibitively slow. 



2. The search 

We limit our search to the simplest class of small- 
universe models, the ones where space is flat and toroidal. 
For this simple case, the matched circles come in diamet- 
rically opposed pairs, specified by three parameters: the 
position, given by (I, b) for the center of a circle, and its 
angular radius a. A much more ambitious six-parameter 
search corresponding to the general case will be presented 
by Spergel et al. [37]. 

To search for matched circles, we define a family of 
difference maps. For a given axis rij, the map is defined 
through 

d^^o-^W (io) 

where Sy denotes the reflection of rij in the plane whose 
normal is rij (i.e., is given by equation (5)). In other 
words, the map D is just the original WMAP map minus 
its reflection about the plane normal to rij. Note that 
equation (10) and equation (4) are closely related: So 
is simply the smallest average of a squared D-map. In 
the ideal case of a perfectly matched circle pair, it would 
manifest itself as a pair of D = circles in a D-map that 
otherwise fluctuates with about twice the variance of the 
original WMAP map. Note that in addition, there will 
always be a great (a = 90°) circle with D = lying in 
the reflection plane: here the WMAP map is of course 
equal to its reflection, so the D-map will automatically 
vanish. 

To distill out this matched circle information, we re- 
place each D-map by a single curve d(a; rij) giving its rms 
along circles making a constant angle a with the reflec- 
tion axis, i.e., with rij • ri, — cos a. An example of such 
a curve (for a fixed rij) is shown in Figure 4. In practice, 
the HEALpix pixelization subdivides the sky into 12n 2 
pixels, where the n is a power of two determining the reso- 
lution. We compute d(a) in n angular bins equispaced in 
cos a, which corresponds to a bin width about 0.87 times 
the pixel side (47r/12) 1/2 /n for a = 30°. The advantage 
of this binning scheme is that each bin gets contribu- 
tions from about the same number of pixels, 12n. We 
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0.2 



0.4 0.6 

cos a 



FIG. 4: An example of the curve d(a; n;) that we used to search for 
matched circles. This case is for the reflection axis correspond- 
ing to (l,b)= (—110°, 60°). A pair of perfectly matched circles of 
angular radius a would give d{a) = 0. The drop towards zero on 
the left hand side is caused by a great (a = 90°) circle being its 
own reflection, whereas the high values around a ~ 15° are caused 
by residual foreground contamination. 



performed our analysis with n — 256, which took a cou- 
ple of weeks on a 2 GHz linux workstation. For a given 
potential circle center n^, the curve d(a; n^) will thus os- 
cillate randomly as a function of a with roughly constant 
variance, as shown in Figure 4. Apart from pixelization 
effects, this curve will drop to zero at any a correspond- 
ing to a perfectly matched circle. Figure 4 shows such 
a deep minimum only for a = 90°, the above-mentioned 
case of the 90° circle which of course equals its own re- 
flection. 

To analyze the results of our search, we compute a 
summary map c? m ; n , defined as 



d n 



,(fii) = min d(a;ni). 



(11) 



If matched circles are present, the summary map should 
show evidence f or these circles in the radius range from 
cto to < a\, and each pixel in this map should correspond 
to the evidence for circles centered on that pixel. 

Figure 5 shows the summary map for circle radii in 
the range 0° < a < 15°. The color scale is seen to 
range over about a factor of two in rf-value (the units 
are /iK). A matched circle would correspond to a pair 
of diametrically opposite pixels with values near zero. 
Not only are no such pixels seen, but the map shows 
that the region around (l,b)= (—110°, 60°) that we are 
particularly interested in given the quadrupole-octopole 
anomalies is in no way unusual compared to the rest of 
the map. The larger values seen near the equator of the 



map are caused by residual foreground contamination, 
which from the definition of the d m ; n maps propagates 
about one circle radius from the Galactic plane. Whereas 
this foreground contamination can mask true matched 
circles centered in this region, it can clearly never create 
false positives, since it increases rather than decreases 
the (2 m i n -value. 

We computed analogous maps for four other ranges of 
circle sizes a, 15° - 30°, 30° - 45°, 45° - 65°, 65° - 85°, 
again finding no evidence of matched circles. For com- 
pleteness, we also computed a final map for the a-range 
85° — 90° obtaining roughly zero, since all 90° circles are 
matched with themselves. 

In conclusion, we find no evidence for the matched cir- 
cles that arc predicted for a simple toroidal universe. 



IV. CONCLUSIONS 

We have investigated anomalies reported in the Cosmic 
Microwave Background maps from the WMAP satellite 
on very large angular scales. There are three independent 
anomalies involving the quadrupole and octopole: 

1. The cosmic quadrupole on its own is anomalous 
at the l-in-20 level by being low (the cut-sky 
quadrupole measured by the WMAP team is more 
strikingly low, apparently due to a coincidence in 
the orientation of our Galaxy of no cosmological 
significance) . 

2. The cosmic octopole on its own is anomalous at the 
l-in-20 level by being very planar. 

3. The alignment between the quadrupole and oc- 
topole is anomalous at the l-in-60 level. 

Although the a priori chance of all three occurring is 
1 in 24000, the multitude of alternative anomalies one 
could have looked for dilutes the significance of such a 
posteriori statistics. 

However, in the context of constraining cosmic topol- 
ogy, these three effects are not merely three arbitrary 
ones among many. Indeed, the simplest small Universe 
model where the Universe has toroidal topology with one 
small dimension of order half the horizon scale, in the di- 
rection towards Virgo, could explain 1, 2 and 3. In order 
to test this hypothesis, we applied the S'-statistic and the 
circle test on the WMAP data, ruling out this model. In 
other words, we have ruled out the "plain bagel" small 
universe model. 

Our results also rule out other models that predict 
back-to-back matched circles. However, they do not rule 
out the recently proposed dodecahedron model of [35]: 
although this model predicts six pairs of diametrically 
opposed circles of radius about 35° , the circles have a 36° 
twist relative to their twin images, thereby eluding our 
search method. After the original version of this paper 
had been submitted, a more thorough analysis by Cor- 
nish and collaborators [36] confirmed our findings and 
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FIG. 5: Result of search for matched circles of radius 0° < a < 15°. A perfectly matched circle would show up as a pixel with zero 
temperature at the position of the circle center, whereas the map above shows no pixels below 83/iK. Note that this map is parity-symmetric, 
i.e., the temperature at n equals that at — n, although this symmetry is obscured by the Aitoff projection used. 



improved them to rule out this and other twisted back- 
to-back models as well. 

A maximally ambitious six-parameter "everything 
bagel" circle search, corresponding to the general case 
of arbitrary topologies, is currently being carried out 
by Spergel and collaborators, and will be presented in 
a forthcoming paper [37]. This should provide decisive 
evidence either for or against the small universe hypoth- 
esis. If this circle search confirms our finding that small 
universes cannot explain the anomalies, we will be forced 
to either dismiss the anomalies as a statistical fluke or to 
search for explanations elsewhere, such as modified in- 
flation models [21-26]. Even the fluke hypothesis might 
ultimately be testable, since it may be possible to im- 
prove the signal-to-noise of the large scale power spec- 
trum beyond the WMAP cosmic variance limit by em- 
ploying cluster polarization [38, 39] or weak gravitational 
lensing [40] techniques. 
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APPENDIX A: ROTATION 

In this appendix, we review the rotation properties of 
spherical harmonics that were used in Section II B and 
Section II C. This material is well-known. However, since 
we found it rather time-consuming to assemble these ex- 
plicit results from the literature (which contains a variety 
of notational conventions and is mainly geared towards 
the general quantum case where functions are complex 
rather than real), we summarize the results here for the 
benefit of the reader interested in performing similar cal- 
culations. Fortran77 software implementing this is avail- 
able from the authors upon request. 



computed in equation (1). The axis n mentioned in 
Section II B thus does not involve cj> and is defined by 
n = (sin cos a, sin sin a, cos 9). The y-rotation multi- 
plies the vector of (2£ + 1) a£ m -coefficicnts for the £ th 
multipole by the Wigner matrix ~D l (6) given by [43-45] 



D 



AO) 



m'm{£-\-m ,£-\-n} 

£ 

k — m ax { , rn + n } 



2k—m—n / Q\2£-\-m-\-n—2k 

I sin - ) 

2) \ 2) 

y/{£ + m)\{l - m)\{£ + n)\{£ - n)\ 



k\{£ + m - k)\{£ + n- k)\{k -m-n)V 
Transforming to our real-valued spherical harmonic ba- 
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1. Real- valued spherical harmonics 

When a function ip(n) is real- valued, the corresponding 
spherical harmonic coefficients obey ai- m = (— l) m a| m . 
The need to work with complex numbers is conveniently 
eliminated by working with real-valued spherical har- 
monics, which are obtained from the standard spheri- 
cal harmonics by replacing e vm ^ in their definition by 
V2sinm^, 1, V2 cosm(j> for m < 0, m = 0, to > 0, re- 
spectively. In other words, the vector of 2£ + 1 complex 
numbers ag m (m = —£, ...,£) specifying the £ th harmonic 
is replaced by a vector of 2£ + 1 real numbers at m given 
by \[2 Im ci£ m , a o, V2Reae m . This mapping from a£ m 
to de m corresponds to multiplication by a unitary matrix 
U; for example, the £ = 3 case corresponds to 
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By unitarity, J2e \ a em\ 2 — ^2g0^ m , and it follows that 
the desired right-hand-side of equation (1) is simply 
The real-valued spherical harmonic coefh- 
observed by WMAP are listed in Table III for 
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for the quadrupole case, and the corresponding matrix 
for rotations around the z-axis is 
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For the octopole case, equation (A3) is replaced by the 
block-diagonal matrix 



where 
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2. Rotations 

An arbitrary three-dimensional rotation is specified by 
three Euler angles (4>,9, a), defining a rotation by a 
around the z-axis followed by a rotation by 9 around 
the y-axis followed by a rotation by <f> around the z- 
axis. The z-rotation by <fi has a trivial effect on the 
spherical harmonic coefficients, simply multiplying ag m 
by a phase factor e*" 1 ^, so the Euler angle <j> will not 
affect the absolute value \ai m \ and hence the quantity 
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The generalization of equation (A4) to arbitrary mul- 
tipolcs is obvious. 



